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Abstract 

We explore constraints on inflationary models employing data on large scale 

structure mainly from COBE temperature anisotropies and IRAS selected 

galaxy surveys, taking care not to apply linear perturbation theory to data in 

the non-linear regime. In models where the tensor contribution to the COBE 

signal is negligible, we find the spectral index of density fluctuations n must 

exceed 0.7. Furthermore the COBE signal cannot be dominated by the tensor 

component, implying n > 0.85 in such models. The data favors cold plus hot 

dark matter models with n close to unity and ^hdm ~ 0.20 — 0.35. We 

present realistic grand unified theories, including supersymmetric versions, 

which produce inflation with these properties. 
98.65 Dx,98.80 Cq,12.10 Dm 
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I. INTRODUCTION 



The inflationary universe scenario |T|] provides an attractive resolution of the well known 
"horizon", "flatness", and "causality" puzzles encountered in the standard big bang cos- 
mology. This scenario is most easily realized within the framework of both ordinary and 
supersymmetric grand unified theories (GUTS). In the simplest realizations of inflation, the 
density fluctuations are Gaussian and close to the Harrison- Zeldovich form, and the back- 
ground density equals the critical density pc- Primordial nucleosynthesis implies that less 
than 10% of the background density is composed of baryons, and so the bulk of matter is 
"dark" , presumably in the form of relic elementary particles. 

Almost a decade ago it was pointed out by Shafi and Stecker that two types of 
dark matter, cold and hot, would provide, within the inflationary context, an elegant way 
of understanding observations which indicated a surprising amount of clustering on larger 
than cluster scales. Since only the cold component clustered on smaller scales, and the hot 
component only clusters on larger scales, such a universe would have an enhanced large 
scale clustering power. Examples of GUTS containing cold plus hot dark matter (C+HDM) 
were also presented. The implications of this picture for microwave background anisotropies 
were worked out in 1989 0,0. Mass functions bulk streaming motions 0,^, cluster 
number densities and correlation functions, as well as "great attractors" were also found to 
be compatible with data. 

Recent data on large scale structure from a variety of sources, particularly COBE and 
IRAS selected galaxy surveys, provide additional strong support to this remarkably "simple" 
scenario for structure formation. The COBE group |§| found an anisotropy amplitude which 
they characterized by an averaged quadrupole moment amplitude 5T/T = Qrms-ps /T = 
6.2 ±1.5 X 10^^. As we have emphasized 0, the earlier predictions for C+HDM models 
(3/4 CDM, 1/4 HDM) with negligible baryonic matter content were 7.8 x lO^^/bg, while for 
CDM alone they were significantly smaller (4.7 x 10^^/6), where fog =(rms mass fluctuation 
on the scale 8 /i~^Mpc)~^. Including a 5% baryon density (implied by primordial nucleosyn- 



2 



thesis) modifies these predictions to 8.2 x 10~^/bs for this C+HDM mixture and 5.1 x 10~^/bs 
for CDM. Observationally derived values of bg seem to fall in the range 1.3 < b^ < 2.5, mean- 
ing the COBE result was in remarkable agreement with C+HDM models. COBE also lent 
some weak support for the inflationary power spectrum; the data are consistent with the 
P{k) oc /c" spectrum with n = 1.1±0.5 Following the COBE result other pieces of 
evidence pointing to the C+HDM model were discovered: compatibility with early quasar 
formation [0,^, compatibility with galactic correlations and velocities, |]12|-Q, APM cor- 



relations [fl4| , p!5| quiet local Hubble flow "cosmic mach number" ^7j, and counts in 



cells Thus, there seems no doubt that the C+HDM model is a candidate worthy of 

serious consideration. 

While we know that the primordial density power spectrum P{k) oc k with C+HDM 
provides a good fit to the data, it has long been known that the inflationary scenario does not 
quite yield this simple form. In recent studies, these correction factors have been exploited 
to yield power spectra P{k) oc k"' with n significantly less than 1 as a means to repair the 
relative large/small scale problems of the pure CDM model. In this paper, one of our tasks 
will be to undertake a study of the range of allowed n. [Note that in this work we will 
only consider n<l, although inflationary models where n exceeds unity can be constructed.] 
Since we do not precisely know the HDM fraction of the universe, we will explore the two 
dimensional parameter space n — ^hdmi and examine the constraints imposed by the data. 

We do this by comparing the model predictions with data from COBE, the power spec- 
trum of IRAS galaxies, IRAS counts-in cells, and bulk velocities from the POTENT analysis. 
We will also consider requirements for the early formation of quasars. After fitting the nor- 
malization and the bias factor for IRAS galaxies, we calculate foi' each model. We then 
compare the results of the models and present contour plots in the n — ^Ihdm fraction 
plane. We also consider the possibility that some of the COBE signal was produced by the 
long wavelength gravity waves generated during inflation. 

Because of uncertainties in the simulation of structure formation in the non-linear gravity 
regime, we will concentrate chiefly on structure which is still described by linear perturbation 



theory (scales > 20h~^ Mpc). This avoids the difficulties inherent in identifying galaxies and 
clusters and separating out complicated dynamic effects, e.g., the velocity bias of galaxies. 
Our main aim is to explore a large region of parameter space, and try to identify that part 
of it for which running more detailed computer simulations would make the most sense. 

Finally, if we are to take seriously the whole picture of inflation and C+HDM we need 
to identify plausible models of inflation which are compatible with our analysis of the large 
scale data. Here we will consider only "realistic" models of inflation. By this, we mean 
models in which the inflating field is a part of a larger theory which contains the following 
elements: 

• neutrinos with masses in the eV range. 

• a cold dark matter candidate. 

• a high energy particle physics structure which is compatible with low energy physics 
and its hints about higher symmetries. 

• a successful baryogenesis following inflation. 

Without these elements we find the infiation model to be somewhat ad hoc. 

The plan of the paper is as follows. In the next section, we will describe our testing of 
infiationary models with large scale structure data in the linear perturbation regime using 
a analysis. This will turn out to strongly limit the power spectrum. In section III we 
will consider some general constraints from data which come from structure in the non- 
linear regime. This will limit mainly the HDM fraction. In section IV we consider models of 
infiation which satisfy all of our requirements for "realistic" infiation models. In particular we 
present some new examples of chaotic infiation based on supersymmetric grand unification. 
We end with some general conclusions in section V. We have attempted to make this paper 
somewhat self-contained for a more general audience. 



4 



II. CONSTRAINTS FROM LARGE SCALE STRUCTURE DATA IN THE 

LINEAR REGIME. 



In this section, we will describe data and our analysis of theoretical predictions for 
these observations, but first we will discuss a few general issues concerning the models and 
the testing strategy employed here. 

First of all, the problem with drawing strong conclusions about any given model of 
structure formation is that there are so many parameters to vary that we have a multi- 
dimensional parameter space to explore. We will choose to take best guess values of three 
of them, the baryonic fraction flbaryon = Pbaryon/Pc, the cosmological constant A, and the 
Hubble constant, Hq = lOO/i km sec"^, with h observationally constrained to the values 
0A<h<1.0. Constraints on the age of the universe from globular star clusters, nuclear 
cosmo chronology, and white dwarf ages imply that high values of the Hubble constant are 
forbidden, i.e. h < 0.6 in an f2 = 1 universe. Thus these models will be allowed only if 
the observations eventually settle into a more restricted range 0.4 < h < 0.6. Since many 
quantities vary as h"^ there is still some freedom despite this narrow range of Hubble constant. 
In the present analysis we will use only the central value, h = 0.5. 

Primordial nucleosynthesis strongly constrains the baryon density as 0.010 < 
flbaryonh"^ < 0.015. Usiug h = 0.5, wc cau express this 95% confidence constraint as 

nbaryon = 0.05 ±0.01 (2.1) 

We note that allowing for the uncertainty in h, the baryonic fraction could range from 
^baryon ~ 0.03 — 0.9. We intend to explore this uncertainty in the Hubble constant and 
baryon fraction and its implications in a future publication. Once we have fixed the baryon 
density, the other densities can be completely specified by the hot dark matter density fin dm 
as 

^CDM = 1 — ^baryon — ^HDM- 

The hot dark matter fraction (combined with the bubble constant) also specifies the neutrino 



mass. If we have one flavor of neutrino whose mass is in the eV range, usually taken to be 
Ut-, then 



.23 eVy \ h J 

We will set the cosmological constant (A) to zero. There is evidence to support this 



choice. The local (within 60 h^^ Mpc) velocity fleld implies values of Q close to unity pi 
Furthermore, in a A dominated universe, there seems to be too few gravitational quasar 
lensing events, and the bulk streaming velocities are too small. 

The growth of density fluctuations is affected by the dynamics of the matter content, 
producing a scale dependent modiflcation of the density fluctuation power spectrum. The 
relative growth as a function of scale is discussed in Appendix A, and the results are sum- 
marized in Figure 1, where we present the transfer functions in Fourier space as a function 
of Fourier wavenumber k = 2tt/X. In flgure 1, we see that increasing the neutrino fraction 
decreases the amount of growth, and hence the amplitude, on small scales (large k). We also 
note that there are some modiflcations of the shape of the spectrum on quite large scales. 
For example the transfer function with ^hdm = 0.3 is slightly steeper at A; ~ 0.07/i Mpc~^ 
than that of either VLhdm = 0.0 (CDM) or fi^mf = 0.5. 

In order to do our comparison in the most unambiguous way, we will conflne our attention 
mostly to the regime where linear theory is appropriate. In the past, it was common to 
normalize linear power spectra by use of the fact that the rms optically selected galaxy 
density fluctuations are 6N/N = 1 on a scale of 8 h^^ Mpc. In the most naive case, where 
one assumes that galaxies (light) trace the mass distribution, one would set the rms mass 
fluctuation also equal to one on this scale. More recently, however, it has become apparent 
that the galaxies are more strongly clustered than the mass, i.e., they are biased tracers 
of the mass. The usual method for dealing with this complication is to assume that there 
is a linear relation between the optical galactic density and mass fluctuations using a bias 
parameter b, which is independent of scale b = {6N/N)/{6M/M). The mass fluctuation 
cr(-R) in a sphere of radius R is calculated in linear theory by 
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where ji{x) is the first spherical Bessel function and the term in brackets is the Fourier 
transform of a sharp edged sphere of radius R. Clearly, when a perturbation amplitude 
approaches unity, linear theory is no longer appropriate. Even with b as large as 2, this 
implies a{8 h^^ Mpc) = 0.5, which is still quite non-linear. Using the spherical collapse 
approximation (see Appendix B) we estimate that for a = 0.5, we are already highly con- 
taminated by non-linear effects. If we consider scales for which (t{R) < 0.4 we estimate that 
the non-linear corrections will be <30%. 

A look at figure 1 shows that the greatest difference between the models considered here 
occurs on the small scales. Limiting our range to k < 0.3h Mpc~^ implies that our testing 
will have a somewhat weakened ability to discriminate between models with different HDM 
fractions. In section III we will consider some constraints from non-linear structures to help 
us pin down the dark matter fraction. 

We will first give a brief description of the large scale structure data, followed by a 
description of our calculations for different models. 



Here we will discuss the particulars of the data we are using. We explain our reasons for 
choosing the data and method of interpretation of this data. 

1. COBE data 

The large amplitude of the COBE measured temperature fluctuations is characterized 
by the extrapolated quadrupole moment Qrms-ps- This amplitude was a factor of 2-3 larger 
than predicted in the usual CDM models fit to galactic structure, and so gave strong support 
for the C+HDM models |^. However, various authors seem not content to use the COBE 
analysis of the amplitude. Instead they choose to use the sky variance at 10° which, when 
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A. Large Scale Structure Data 



used with a Gaussian beam, implies a somewhat smaller amplitude for the fluctuations (e. 
g., ||r5| , |22[] corresponding to Qrms-ps = 15.3yuK). The COBE beam pattern, however, is 
only approximately a Gaussian shape. More careful analysis of the COBE results I^S] using 
the actual beam pattern seems to confirm the original higher value of the amplitude. In 
addition, use of the correct beam pattern reduces the variance of the fit amplitude, so the 
best fit n = 1 COBE amplitude now corresponds to Qrms-PS = 17.1 ±2.9 /iK. Although it is 
customary to quote the results in terms of the quadrupole moment, this corresponds to only 
the very small wavenumber end of the spectrum. The best fit amplitude of the quadrupole 
moment will be somewhat dependent on the value of n used in the analysis. In order to 
find a better quantity than the quadrupole, Wright et al. recommend normalizing the 



amplitude to the hexadecupole ST^^ = 12.8 ± 2.3 /iK, when n 7^ 1, as in the infiationary 
models we consider here. This yields a best fit amplitude less dependent on the value of n. 
In passing, we note that the COBE results should be taken seriously, as a balloon exper- 



iment from the MIT/GSFC/Princeton collaboration |2^ sees the same temperature corre- 
lation function as COBE and find similar values for the fiuctuation amplitude, which they 
specify with Qrms-PS- However, they currently only have data from the northern hemisphere 
and their limits on Qrms-PS are not as constraining as those for COBE. In the future they 
plan to cover the southern hemisphere as well, which would improve the limits of the COBE 
power spectrum exponent, which are currently n = 1.1 ± 0.5. 

Ideally, we would like to include the detections of smaller scale anisotropy experiments in 
this analysis. There has been a wave of new detections of temperature anisotropics on 1° — 2° 
scales. It is not clear that all of the detections are giving a totally consistent picture, and 
even different scans with the same instrument give different results. A possible explanation 
is that the systematic errors these experiments face are quite complicated. We expect 
that these uncertainties will be resolved and that we will know the amount of degree scale 
anisotropy which exists with some precision. We also point out that C+HDM models, when 
normalized to COBE produce the same degree scale anisotropics as similarly normalized 
CDM models. The only differences occur for anisotropy measurements <^ 1°. Since the 
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degree scale anisotropies are all in the right ballpark for Q = 1, n ^ 1 models normalized to 
COBE, we take this as an encouraging sign for the models we consider in this paper. 



2. The IRAS Large Scale Survey of Galaxies 

The IRAS survey of galaxies done by the QDOT (Queen Mary-Durham-Oxford- Toronto) 
collaboration extends as deep as several hundred Mpc, approaching the smallest scales 



observable by the COBE satellite. Combining the COBE data with the IRAS survey of 
galaxies thus covers the whole range of scales where the fluctuations can be described by the 
linear theory. The QDOT IRAS survey has measured redshifts for 1 - in - 6 IRAS galactic 
sources (1824 galaxies) with IRAS fluxes > 0.6 Jy. The IRAS selected galaxies seem to be 
more uniformly distributed than optically selected galaxies and may therefore give a fairer 
representation of the universe. By concentrating on measuring only 1 in 6 galaxies, the 
QDOT collaboration obtained a deep sparse sample out to a depth larger than that of the 



Berkeley IRAS survey which measured redshifts for all IRAS sources above a flux of 1.2 
Jy. Using the redshift of the source as a distance indicator, combined with angular position 
data, one has a three dimensional picture of the galaxy distribution. This distribution can be 



analyzed to directly extract the power spectrum of density fluctuations The results 

of the Feldman, et al. |2^ analysis is shown as the power spectrum data in Figure 2. 

We will use this QDOT IRAS data to test models with theoretical power spectra given 

by 

Pthik) = Ak-[T{k)]\ (2.4) 

For n ~ 1 and the transfer functions T{k) in figure 1, we can see that the power spectra 
go like P{k) ~ on very large scales and like P{k) ~ k~^ — k~^ on very small scales, with 
a peak somewhere around k ~ 0.02 — 0.10. In order to get the proper normalization for 
the amount of structure on scales up to the power spectrum peak scale, one must measure 
the power on scales larger than the power spectrum peak scale. In the mass fluctuation 
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integral (equation 2.3) there are significant contributions to the mass fluctuations in 20 Mpc 
spheres coming from 100 Mpc scales {k ~ 0.03/i/Mpc). Thus we feel that to test power 
spectra of the type considered here, we require a survey out to a depth at least as large as 
the QDOT survey. A comparison of the QDOT and 1.2 Jy Berkeley power spectra, (see Ref 
P7|), suggests that the Berkeley survey may still be too small to be seeing all of the large 
scale power. 

On large scales, the power spectrum is a better measure of clustering power (i.e., has 
smaller errors) than the galaxy correlation function |2^. However, for scales < 100h~^ Mpc, 



the power spectrum is not as good an indicator. Since the largest differences between the 
models under consideration here occurs on smaller scales, we would like to supplement our 
analysis with a different estimator of power on small scales. To this end we use the "counts 



-in - cells" statistic from Efstathiou, et al., Ref. ||2^. Here the redshift space of IRAS 
galaxies is cut up into near cubical cells of side length i and the number of galaxies in each 
cell is counted. The cell to cell variance of galactic number is a direct indication of the 
underlying density fluctuations of length ~ £ in which the galaxies reside. The statistic is 
usually denoted by the symbol <j'^{i)- The values given in Ref. |]25| are cr^(^) = 0.42 ± 0.07, 



0.26 ±0.05, 0.21 ±0.05, and 0.047 ±0.024, for £ = 20, 30, 40, and 60 h'^ Mpc, respectively. 
They also give a value for i = lOh^^ Mpc, but this represents a fluctuation which is strongly 
contaminated by non-linear effects, and hence not appropriate for our analysis. The QDOT 
counts in cells data is presented in flgure 3. 

3. POTENT bulk flow velocities 



We also use the bulk flow velocities from the POTENT analysis |2^. They represent 
the the rms velocities of spherical regions of radius R. The velocity has flrst been filtered 
with a Gaussian of width R = 12h^^ Mpc to reduce noise. This data is shown in figure 
4. The POTENT analysis has been done using the IRAS galaxies of the Berkeley survey, 
and so suffers from the problem that enough of the universe has not been sampled to get a 
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fair estimate of the velocities. This is exacerbated by the fact that velocities are even more 
sensitive to the very large scale power than the mass fluctuation. However, we would still 
like to find a way to use these velocities because they do not depend on the bias. (Velocities 
are generated by the gravity of density fiuctuations.) Since we have velocities from one local 
patch we will include the cosmic variance of the predicted velocity for any given patch of 
the universe in our analysis. Our treatment will be described in the next section. 



B. The Test 

The reduced statistics is calculated by comparing the set of N predictions yf^ and 
observations yf"^ according to the following formula 

x' = j/-T{^!^-J-) . (2.5) 

where a°^'^ are the standard observational errors. Nct.o.f. equals the number of observations 
N minus the number of fitted theoretical parameters in the model. In our analysis we use 
the 36 values of the IRAS P{k), the 4 IRAS counts-in-cells, the 5 POTENT velocity values, 
the bj determination, and the COBE amplitude, which we will count as two points. In the 
analysis of n = 1 models in ref. |jl8|, the COBE data is counted as two points, one for 



the sky noise at 10° and another for the value of Qrms-ps- Because of the complications 
of the non-gaussian beam pattern we have avoided this procedure, but we still choose to 
weight the hexadecupole amplitude as if it were two points in the analysis. This forces the 
amplitude of density perturbations to be slightly closer to the COBE normalized amplitude 
than by weighting it as only one point. The results of the analysis is quite similar if we were 
to weight COBE as only one data point. Thus we have a total of = 48 data points which 
contribute to the sum. 

For each value of n and VLhdm we find the values of and hj which minimize x^(^8; 
for each model. This is effectively a least squares fitting of and 6/ for each model, so the 
number of degrees of freedom is 2 less than the number of points, Nd.o.f. = 46 in our 
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formula. The confidence levels of rejecting the hypothesis that the data fits the model are 
found by integrating the normalized distribution. The probability of getting > 1.38 
is < 5%, so we can reject such models with 95% confidence. We will now proceed to discuss 
how we compute the yj^. 



C. Model Predictions From Linear Theory 

Here we will explain how we compare our theoretical power spectra Pth{k) to the data 
presented in the previous section. First of all, we can calculate the hexadecupole (5T4 (mea- 



sured by COBE) with the formula given in Ref. for the coefficient of the fourth spherical 
harmonic for each of the power spectra we are considering. However, we can find an effec- 
tive wavenumber k^ff for which the amplitude of P{k^ff) is directly proportional to the 
amplitude of the hexadecupole. We have found that the value fee// = 1-05 x 10^'^/i/Mpc 
accurately characterizes the hexadecupole moment for the limited range of n we are studying 
(0.70 <n< 1.00). 

The COBE experiment cannot distinguish between temperature fluctuations generated 
by the Sachs- Wolfe effect of scalar (density) fluctuations and tensor (gravity wave) fluctua- 
tions, both of which are products of inflation. This signal confusion is worsened by the fact 
that the ratio of the moments generated (at least on COBE scales) for scalar and tensor 
contributions are nearly independent of scale pTj- p^j . The overall amplitudes of the mo- 
ments for scalar and tensor modes, however, can be quite different, and one must construct 
a speciflc model of inflation and then evaluate the density perturbation and gravity wave 
amplitudes during inflation. In some models of inflation, such as "new" inflation [||], the 
contribution of gravity waves to the COBE signal is negligible, and the COBE signal relates 
only to the scalar density perturbations. 

Models of inflation which produce a signiflcant amount of gravity waves cannot be sum- 
marized with a universal formula. However, there is a relation derived from the toy model 
"power law inflation" [Bl] which relates the gravity wave (tensor) contribution to the COBE 
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signal {AT/T)t to the contribution from the density (scalar) fluctuations signal {AT/T)s 
via the spectral exponent of the density power spectrum n p2| , p2| , |35| : 

This relation also is reasonably accurate for chaotic inflation models. Since the COBE 
signal is a quadrature sum of the multipole moments, decreasing n decreases the fraction of 
COBE signal due to density waves, and thus implies a smaller amplitude for Pth{k). Thus 
the amplitude of the COBE quadrupole moment is — 7n times larger than would be 
expected from the density component alone. We will consider models with gravity waves 
produced according to this formula. 

To simulate the IRAS power spectrum Pi{k) we have to apply a couple of correction 
factors to our theoretical Pth{k). First of all the distribution of galaxies in redshift space 
on large scales appears more clustered because of the doppler contribution to the redshift 
from velocity perturbations [Q. With a biased galaxy distribution (bj) this correction can 
be made by multiplying with the factor 



Pik) 



2 1 



Pik) (2.7) 



However, on smaller scales there is the opposite effect; the doppler shifts from the peculiar 
velocities of galaxies become large in comparison to the Hubble velocities and in fact wash 
out this redshift clustering effect. This effect can be described with a velocity dependent 
factor [^] 



P(k) ^ ^ (2.8) 



''V 



where = 4.4 Mpc for IRAS galaxies. 

Applying these corrections to a linear power spectrum seems to give reasonable agreement 
with the power spectrum of "galaxies" in N-body simulations. For example, Feldman, et al. 
P7| present such a power spectrum (from ref. |TB|) for a model with 30% HDM. The N-body 
spectrum becomes slightly higher than our linear spectrum for A;>0.15 presumably due to 
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non-linear corrections. The worst disagreement (<30 %), occurs at the small wavelength 
end k = 0.2h Mpc~^ of the QDOT IRAS P{k). This error is still much smaller than the 
QDOT error bars, so we believe our procedure is relatively insensitive to non-linear effects. 

To compare our predictions with the counts in cells data we calculate the mass variance 
in a spherical volume of radius R = (3£^/47r)^/^, i.e., the radius which encloses the same 
volume as i^. We have checked that this procedure gives the same results as a cubic cell of 
length i by computing the variance in a cube and a sphere of equal volume using a CDM, 
n = 1 power spectrum. The results were different by only a few percent. 

We calculate the counts in cells directly for each model, using the same redshift space cor- 



rection as appropriate for the IRAS power spectrum. Peacock has suggested a different 



method for comparing counts-in-cells data to P{k). He noted that one can write 

2nWil) = P{k,ff)klff (2.9) 

where the value of fcg// depends on the spectrum. This is the method used by Taylor and 
Rowan- Robinson |]TB[ for their analysis. It is instructive to calculate /cg// for a typical 



model just to illustrate the wavenumbers being probed with the counts-in-cells measure- 
ments. For a n = 1, fijy = 0.25 model, the values for k^ff are given in Table 1. 

Thus we see that for accurate count-in-cells numbers for these volumes, one needs to 
accurately get the power spectrum out to A;<0.02. 

Note we are not considering the datum from cells of length I = lOh^^ Mpc, as this point 
is strongly in the non-linear regime. In appendix B we estimate that even for i = 20/i^^ 
Mpc, the non-linear effects may be as large as ~ 30 %, which is roughly the same size as 
the disagreement between our linear P{k) and the N-body simulation P{k) at the smallest 
wavelength we consider, k = 0.2/;, /Mpc. For this reason we do not consider estimates of the 
power on scales much below 20h^^ Mpc. (The largest wavenumber from the IRAS QDOT 
P{k) is kmax = 0.195 h Mpc^^) 

To do the testing, we first must calculate the linear theory power spectra for all models 
with a primordial power spectral index 0.70 < n < 1.00 and < < 0.5. We have 
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limited n to be < 1 because we are considering only grand unified models of "new" and 
"chaotic" inflation. We calculate the power spectra for these models in steps of 5n = 0.02 
and = 0.05 using our form for P{k): 

Pth{k) = Ak^[T{k)]\ (2.10) 

We will vary the amplitude A by a factor of 2^^ away from the COBE implied central value in 
20 logarithmically spaced steps. However, since the definition of A is author dependent, we 
will discuss our results in terms of the parameter feg for each model. Thus we are effectively 
varying fog in a range centered on the COBE best fit fog- 

The value of 6/ is somewhat more constrained than values of 6g. For the IRAS galaxies, 
separate determinations of 6/ have been made by comparing their velocities and distributions. 
These dynamical tests yield the 95% confidence values values hj = 1.16 ± 0.42 ( Ref. p8| ) 
and bi = 1.23 ± 0.46 (ref. |9|). However, the POTENT analysis ^ of the Berkeley IRAS 



galaxies finds that the 95% confidence interval for bj = 0.5 — 1.3. We therefore combine 
these measurements to say that bj = 1.1 ± 0.3 to cover the 95% confidence overlap region 
of the bj determinations. We consider the measurement of 6/ = 1.1 ± 0.3 to be a bona fide 
data point which we add to our test. As is usual, we will assume a constant bias factor 
independent of scale and we will allow bj to vary in 7 linearly spaced steps between 0.8 and 
1.4. We calculate the for sach of the values of bj and 6g, and find the the set (6/, 6g) 
which gives the lowest value of x^, and plot this minimum in our -- n contour plots. 
Thus we allow each model to put its best face forward in our test. Some of these models 
with bi = 1.00 are shown in figure 5. 

We calculate the POTENT rms velocities by the same procedure as described in ref. 
p!0| . In our comparison of the rms POTENT velocities to theoretical rms velocities, we will 



incorporate the cosmic variance of the velocity field. The idea is that with a Gaussian density 
field, the velocity field will also be Gaussian. The magnitude of the velocity vector will have 
a distribution with 3 degrees of freedom. The variance of the rms velocity is much larger 
than the POTENT velocity errors. For example the 68% theoretical confidence range on 
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the average predicted velocity magnitude (v) corresponds to {{v) — 0.48(v), (v) + 0.32(v)), 
while the POTENT errors are less than ~ 15 %. This is plotted in figure 4 for the n = 0.96, 
25% HDM model. If we normalize an n = 1.00 power spectrum to COBE then the predicted 
velocities are within the POTENT 1 a error bars, regardless of the HDM fraction (see e.g., 
n — 1.00, 50% HDM model in figure 4). As we decrease n below 1.00, we will decrease the 
predicted velocities well below the POTENT values. Thus the upper limit on the predicted 
velocity will be the most relevant quantity for comparing theory to observations. We combine 
this upper hmit in quadrature with the POTENT 1 a error, and use this as a fairer estimate 
of the error bar in our analysis. 

While adding errors in quadrature is strictly correct only for Gaussian errors, the velocity 
distribution is not too far from a Gaussian. We plot several theoretical predictions for (v) 
against the POTENT velocities in figure 4. Even with these huge error bars, we find that 
the predicted velocities for models with significant gravity wave contributions will still have 
trouble matching the POTENT derived velocities. 

D. Results of Linear Data Analysis 

The results of our test are presented in figures 6 and 7 for models without and with 
gravity waves, respectively. Starting with our absolute best fitting model n = 1.00, 30% 
HDM and working outward, we plot 9 concentric curves with confidence levels of .5, 1, 5, 10, 
25, 50, 68, and 95 %. (For reference, the best fit model has formally a probability 4 x 10~^ 
for getting such a good fit by chance - although such small probabilities are meaningless 
for statistical analysis.) When one reaches the level of 1%, areas where the theory does not 
match the data become discernable in the data plots. The first thing to note about the 
graphs is that the overall level of for models with n ~ 1.00 is quite low. This result 
supports earlier claims that n = 1 models normalized to COBE have sufficient power to 
explain "large scale power" apparent in galactic clustering measures. If the measurement 
errors were smaller our test would be a much better discriminator between models. 



16 



We have two scenarios to discuss: models with and without gravity waves. First we 
point out one general trend which is common to all models, and which does not show up 
in the numbers. As we decrease n we decrease the amount of mass clustering power on 
small scales. To increase the galactic clustering power to compensate for this, one needs 
to increase the amount of galactic biasing 67. Since our program finds the best fit 67, we 
will point out that small n models correspond to highly biased models (large 67). We could 
better limit the models if we had more precise information about 67. 

The contours for models with no gravity waves are shown in figure 6. As can be seen 
we can rule out models with at 95% confidence. There is a slight dependence of n on 

the HDM fraction, with the limit being n > 0.67 for models with no HDM and n > 0.72 for 
models with 50% HDM. This is easily understood since models with a lot of "tilt" already 
have little galaxy scale power, and large HDM fractions exaggerate this behavior. The region 
of best fits occurs in a roughly rectangular region 0.85<n<1.00 and 0.1<J1777)m^0.5 That 
these fits are quite good can be seen from the direct comparisons to data shown in figures 
3, 4, and 5 for a few models. 

The contours for models with gravity wave contributions according to equation 2.5 
are shown in figure 7. It is readily apparent that the allowed parameter space is much 
smaller. Here the 95% confidence limit on n rules out models with n<0.85. Again, there 
is a slight dependence of n on the HDM fraction, with the limit being n > 0.83 for models 
with no HDM, and n > 0.87 for models with 50% HDM. The region of best fits occurs in 
a roughly rectangular region 0.94<n<1.00, 0.1<r2777)M^0.5. Including significant amounts 
of gravity waves forces the normalization of the density power spectrum to be much lower, 
which depletes the amount of clustering power. 

We note that for our 95% confidence limit an amount >50 % of the COBE signal is 
attributed to the effect of density fluctuations. This tells us that to fit large scale structure 
one requires that the COBE signal cannot be dominated by inflation generated gravity 
waves. Our best fit models are those in which at least 80% of the COBE signal is due to 
density fluctuations. 
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Our two conclusions from this analysis are that the power spectrum must be close to 
the Harrison-Zeldovich form [n = 1.00) and that the COBE signal must be mostly due to 
density fluctuations. 

III. CONSTRAINTS FROM DATA ON NON-LINEAR STRUCTURES 

As noted in the previous section, we find the best models are those which have 
^i?_D A/^0.1. The statistical confidence of this conclusion is not very high, on the order 
of 10% for a given value of n. This is not surprising, since we have confined ourselves to 
only the largest scales, k<0.2h Mpc^^, where the transfer functions for all the models do 
not differ too much (see figure 1). The best place to discriminate between these models is 
on smaller scales, where we have data only on the non-linear part of the power spectrum. 

A. The Data 

The relation between non-linear structures and the amplitude of the linear power spec- 
trum is quite complicated, and it is non-trivial to extract strong conclusions from this data. 
We will consider two such constraints which we feel can be used relatively safely - constraints 
on the amplitude of cr{8h^^ Mpc) or equivalently bg, and the requirement that quasars form 
early enough to be compatible with observations. The quasar constraint is a lower bound 
on the power spectrum amplitude while a{8h~^ Mpc) is an upper bound (at a somewhat 
larger scale). 

1. High redshift quasars 

The discovery of quasars with high redshifts (about 20 with z > A) was a direct challenge 
to theories of structure formation. One needs a minimum amplitude for density fluctuation 
on galactic scales to account for the quasar population. Efstathiou and Rees (ref. P0[ ) 
considered the formation of quasars in a highly biased (feg = 2.5) n = 1 CDM model. 
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The basic strategy is that in order for a massive black hole to form and power the quasar 
emission, one first requires a host dark matter halo to supply the gravitational potential to 
induce baryonic infall. The number density of structures for a given power spectrum can be 
computed using the Press-Schecter or BBKS techniques. The number density of structures 
depends exponentially on a parameter v given by 

where R is the radius of the initial collapsing region appropriate for the objects in question 
and (Tc is the linear theory value of cr(i?) which corresponds to the gravitational collapse and 
virialization of the object. Using the spherical collapse approximation (see Appendix A) 
(Tc = 3(127r)^/^/20 = 1.69. In this application we use the Gaussian filtered mass fluctuation 

aa{Rf = ^ / dkeP{k)e^'"'^' (3.2) 
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where the subscript "G" is to remind us that we are using a Gaussian filtering function. The 
parameter u then has the physical meaning of the ratio of the overdensity in a collapsed ob- 
ject to the ambient rms overdensities. Thus specifying u tells us the amplitude of the density 
fluctuations on the scale R. Also, since the calculated number density depends exponentially 
on z/, it is hoped that errors in the measured number density will not significantly change 
the value of u. 

Efstathiou and Rees (ref. |^^) estimated that the minimum mass for a quasar halo was 
> 2 X IO^^Mq. They found that a fog = 2.5, n = 1 CDM model could account for the number 
density of quasars at least out to a redshift of z = 5. Haehnelt and Rees (ref. [^) improved 
on this treatment by showing that the b = 2.5 CDM could fit the quasar luminosity function 
at a variety of redshifts. Since both increasing the HDM fraction and reducing n lessen the 
amount of power on small scales, it is obvious to ask whether the amplitude of quasar scale 
fluctuations has decreased below that required for making quasars. Haehneldt (ref. |]rT| ) 



used the techniques of ref. to limit the HDM fraction in n = 1 models and n in CDM 
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models, finding that n > 0.75 in CDM models and the HDM fraction must be ^hdm < -3. 
Schaefer and Shafi (ref. [|T0|) showed that 25% HDM and n = 0.94 — 0.97 is compatible with 
the quasar density out to a redshift of 2; ~ 5. 

The Press-Schecter technique is, however, not without errors. Even if we knew the quasar 
number density perfectly, and hence could deduce the exact value of z/, there would still be 
uncertainties. These uncertainties have been discussed before in a variety of places. Here 



we follow the discussion of reference [^. The uncertainties in u can be found by 



Alnz/ = Alna, - ^^^AlnM (3.3) 
dR dM ^ ^ 

First of all it is not clear what value of o"c to use. While most authors use the value 1.69, 
comparison of Press-Schecter results to N-body simulations imply values of cTc anywhere 
from 1.33 to 1.69 [|^. To bracket these values we can assume = 1.5 ± 0.2. Second 



there are errors in the quoted mass of the object. The theoretical value for the quasar halo 
mass is somewhat uncertain, and the mass of the object derived from the luminosity function 
comes with an additional error. To evaluate the total uncertainty one needs to know the 
value of {da{R)/dR){dR/dM). For the quasar halo mass of 2 x IO^^Mq, the Gaussian filter 
radius is R = OM'^ Mpc. For n = 1 models {da{R) / dR){dR/ dM) ranges from 0.14 for 
CDM models to 0.033 for 50% HDM. Assuming a factor of 3 error in the observational and 
theoretical masses, we then have 

Alnz/ = ±0.13 ±0.04 ±0.04 (3.4) 

for 50%HDM {n = 1) which has the smallest uncertainties of the n = 1 models. We conclude 
that there could be an uncertainty of about 20% in the amplitude of the density fluctuation 
(j{R) derived from the quasar density. 

In order to have the proper quasar density at early redshifts, one needs to have a{O.Qh^^ 
Mpc) > 1.1 today (using (Jc = 1.5). The effect of the hot dark matter on the growth of 
density fluctuations at these scales changes the above threshold value by only a few percent 
|TT[. To forbid models which are unlikely to have quasars form early enough, we make the 
replacement 
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whenever a{0.6h~~^Mpc) < 1.1. We have weighted this contribution to by 20 to strongly 
penahze models without sufficient small scale power. Since there are about 40 data points 
in our analysis an amplitude cr(0.6/i^^Mpc) which is 40% (2 x 20%) smaller than our 
threshold amplitude will cause to be so large that the model will be ruled out by the 
quasar data alone. 



2. a{8h'^ Mpc) 

Since it has become traditional to specify the amplitude of linear theory by bg, some 
attention has been paid to determining the value of bg from observations. There are several 
ways of getting at this quantity which seem to be converging on the range feg = 1.5 — 2.0 
for the = 1 models. We will consider some of the more recent attempts to constrain this 
parameter. 

The easiest way to estimate the density fluctuations is by finding the variance of the 
galactic number density. For the IRAS galaxies we have two estimates of the number density 



fluctuations: one from the QDOT survey [43 



6/(T"''"-'~(8/i-^Mpc) = 0.69 ± 0.09 (3.6) 



which is in perfect agreement with the estimate from the Berkeley survey |44 



6/a"°"-'"(8/i-^Mpc) = 0.69 ± 0.04 (3.7) 

The error on the averaged combined measurement is essentially the same as the Berke- 
ley error, bja"'"'^"'-^'^^"-^ (Sh'^Mpc) = 0.69 ± 0.04. The 95% confidence upper limit is then 
6/cr"°"~'™'^"'"(8/i"^Mpc) < 0.77 If we take an extremely small value for 6/ = 0.5 which is the 
POTENT 95% confidence lower limit, this would imply that a"°"-'^''"''''(8/i-^Mpc) < 1.54. 
If we use the spherical collapse model (see Appendix B) to determine what that means in 
terms of the linear density fluctuation, this implies 
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ct'~ (8/i-^Mpc) < 0.71, or feg > 1-40. (3.8) 

To be extra cautious, we will adopt the constraint 

a'~(8/i-^Mpc) < 0.80, or bs > 1.25 (3.9) 

for a new round of testing. 

A few remarks are in order here concerning the value of fog- Early investigations found 
that CDM models with n = 1 required values of fes ~ 2 — 3 to get agreement with 
galactic velocity dispersion data and correlation functions. However, large scale structure 
data required that CDM have smaller values of bs- It was postulated that dynamic effects. 



especially "velocity bias" might accomodate smaller values of bs- High resolution simu- 



lations 1^ however, find bs > 1.4, despite their confirmation of the existence of the "velocity 



bias" effect. This conclusion seems to also hold for n < 1.00 Ref. also suggests 
that this is true even in models with C+HDM. However, they did not consider the effect of 
the dynamical effects of HDM in their study, and this seems to allow us to use bs < 2 with 



significant amounts of HDM p!2|-p^. No systematic study has been done to determine what 
values of bs are allowed, although bs = 1.5 seems to work with n = 1.00, VLhdm = 0.30. Our 
constraint fog > 1-25 seems to be easily consistent with these case studies. 

Another constraint on the mass fluctuation amplitude comes from cluster properties. 
Since clusters are a few Mpc in size, they are an almost ideal choice for determining feg. 
However, to determine the number density of cluster mass structures analytically, one must 
use a procedure such as the "BBKS" method or the Press-Schecter model. There 



are some uncertainties associated with this technique however, which we have discussed in 
the previous section on quasars. In ref. 0| it was found that 6g = 1.12 — 0.96 for models 
with 0% - 50% HDM, based on i? > Abell cluster number abundance. To arrive at this 
number van Dalen and Schaefer (ref. 0) used the spherical collapse model. The R = Abell 
clusters are the poorest Abell clusters and represent initially weaker density fluctuations. 
The smaller amplitude density fluctuations tend to be highly asymmetric and collapse 
faster than spherical perturbations (see e.g. ref. ||5l|), so a spherical collapse model probably 
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gives bs values which are too small. It was estimated that more accurate values for bg would 
be at least 30% larger, i.e., feg = 1.46 — 1.25, consistent with the adopted restriction 68^1-25. 

Perhaps a more reliable way to study the mass fluctuations in clusters is to select them by 
their X-ray temperatures, as this gives a direct indication of the gravitational mass potential. 
These studies tend to give results consistent with much higher feg, implying fog ~ 2.0 — 2.5 



( 6g = 1.6 - 1.9 [H] and 6g = 2.0, [g. Thus we flnd our restriction fog > 1-25 is, if 
anything, too conservative. 

B. Results of Non-Linear Analysis 

In Figures 8 and 9 we again plot the contours for models with and without signiflcant 
gravity wave temperature anisotropics. The effect of our non-linear constraints is clear. The 
restriction 6g > 1.25 forces the normalization of small fin dm, ~ 1 models to be too low to 
match the large scale structure data. This is symptomatic of the CDM models which, when 
normalized to COBE, have too much small scale power. The effect of enforcing the lower 
limit on the amplitude from quasars is to cut out a triangle of parameter space corresponding 
to large ^hdm and small n. This is symptomatic of ^hdm ~ 1 model problems, there being 
not enough small scale power to explain the early epoch of quasar formation. 

What we are left with is a patch of parameter space which has ~ 1 and ^Ihdm = 
0.3±0.2. This is in agreement with earlier studies and shows that models with VLhdm ~ 1/4 
are significantly better fits to the data than with no HDM. The range of n is roughly the 
same as in the linear analysis, although the non-linear quasar constraint has effectively 
chopped off the low n, high HDM fraction corner of parameter space of the previous best 
fits. Our allowed region of parameter space overlaps the allowed region found by Liddle and 
Lyth, (ref. |]5^). However, their analysis took the non-linear constraints of refs. lll|;E31 



face value so their allowed region is somewhat smaller than ours. They noted however, that 
their allowed region was meant to be suggestive of trends in the data and should not be 
taken literally. On the other hand we are taking pains to be overconservative with non-linear 
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constraints in the hope that our hmits will be firmer. 

Thus we find the following properties which it is desirable for inflation to have. We require 
a density perturbation spectrum which is quite close to a Harrison- Zeldovich spectrum as the 
data do not seem to favor much "tilting" . For models with a negligible gravity wave contri- 
bution to the COBE signal, the best fits occur for 0.9<'ri<1.0. In models with some gravity 
wave contribution, we find an even tighter range of best fit n values, namely 0.94<n<1.00. 
Since we used equation 2.5 to determine this, we see that for n = 0.94 the gravity wave 
contribution to COBE is only 20%. Thus the data favor models for which density perturba- 
tions are >80% responsible for the COBE signal. We also find 0.15<i7//DM^0.45 gives the 
closest fits to the data, which implies we would nominally like one fiavor of neutrino with 
a mass = 3 — 10 eV. (The best fits imply a narrower range ^hdm = 0.20 — 0.35. With 
these attributes in mind we proceed to explore possible models for infiation. 

IV. MODELS OF INFLATION 

Grand unified theories (GUTS) provide the simplest framework for implementing the 
infiationary scenario. Although supersymmetric GUTS are currently more popular, for 
completeness we will also consider the ordinary non-supersymmetric versions. The simplest 
example of the latter is provided by 5*0(10) with an intermediate mass scale. The minimal 
non-supersymmetric SU (5) model is excluded both by the precise determination of sin^ 9w 
and by proton decay experiments. This is just as well from our point of view since, as 
observed a decade ago, |5^, non-SUSY 5*0(10) models with an intermediate (B-L breaking) 
scale Mb-l{^ 10^^ GeV) strongly suggest that the tau neutrino mass is in the eV range. 
The presence of a U{1) axion symmetry not only resolves the strong CP problem but also 
provides the cold dark matter component. 

Two versions of the infiationary scenario, 'new' and 'chaotic', are readily realized in 
GUTS. The spectral index n of density fiuctuations in the simplest realistic models typically 
lies between 0.96 and 0.92, although in some versions of chaotic infiation with SUSY GUTS, 
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n could be as low as 0.88. Values of n much smaller than this are not particularly well 
motivated, both from the point of view of model building as well as observations of the large 
scale structure. 

A. Inflation with Non-Supersymmetric SO(IO) 

For definiteness, let us consider the following breaking: 



50(10) — y 5f/(4)e X SU{2)l x 5f/(2)^ — > SU{3)c x SU{2)l x U{1)y 
Mx Mb-l 



A recent two loop renormalization group calculation involving the gauge couplings gives ||56 
Mx ~ 10^^ - 10^^ GeV and Mb-l ~ 10^^=^^ GeV , consistent with the measured values of 
ac{Mz) and ?,m^dw{Mz). 



A simple version of the see saw mechanism for neutrino masses suggests the hierarchy: 

rriy^ : rriy^ : rriy^ ~ : : O(10^^)m^ (4.1) 

Here the three mass eigenstates z/i, z/2, z/3 primarily consist of z/e, and respectively, and 
we have assumed (see later) that the heavy (right handed Majorana) neutrino associated 
with the third family is a factor 10 (or so) heavier than the other two. 

The MSW interpretation of the solar neutrino data suggests that the z/2 mass is 
~ 10~^'^ — 10^^'^ eV . With mtijnt) ~ 130 — 150 GeV , as suggested by recent analyses of the 
electroweak data, we expect the z/3 (essentially z/,-, with some admixture of z/^ and z/e) mass 
to be ~ 5 — 10 eV . (Without the numerical factor in equation (4.1), this mass would exceed 
the cosmological bound.) Note that according to this simple 5*0(10) example, unless the 

— mixing happens to be tiny, the two neutrino oscillation experiments CHORUS and 
NOMAD should determine whether or not the 'tau' neutrino is cosmologically significant. 

The 5*0(10) model with the above symmetry breaking chain suggests the existence of 
some dark matter in the form of tau neutrino. However, two essential ingredients are still 
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missing, implementation of inflation and a candidate for cold dark matter (CDM). [Recall 
that inflation with only hot dark matter (HDM) does not seem compatible with the observed 
large scale structure, especially galaxy formation.] The simplest way to incorporate CDM 
here is to invoke a U{1) axion symmetry ||59[ broken at a scale around 10^^ — 10^^ GeV. 



Both the axions and neutrinos are then cosmologically significant. It may be useful to 
reiterate how this has come about. In non-supersymmetric SO (10), an intermediate scale 
is needed to bring about consistency with the measured value of sin^ 9{Mz) as well as with 
the lower limits on the proton lifetime. This forces the gauged B — L symmetry to break 
at an intermediate scale which, for the above chain, is about 10^^^^ GeV . Coupled with the 
see saw, this strongly suggests that the tau neutrino mass (~ ml/MB-L) is in the eV range. 
That is, the neutrino is a significant component of the dark matter. Cold dark matter is 
then needed to reconcile the infiationary scenario with observations related to large scale 
structure. 

As far as infiation is concerned, the most straightforward scenario is realized by intro- 
ducing a weakly coupled gauge singlet field a la Shafi and Vilenkin |^. The part of the 
potential which drives (new) infiation is given by 

^('^) = ^<^'M^-^) (4.2) 

where M denotes the vacuum expectation value (vev) of (j). The quantity A can be reliably 
estimated by considering the contribution of the scalar perturbations to the microwave back- 
ground quadrupole anisotropy and identifying it with COBE's determination of Qrms-ps- 
[Note that for the potential in (4.2), the spectral index n ~ 0.94, and the tensor contribution 
to the anisotropy is negligible.] 

One has (the subscript S denote the scalar contribution) 

'AT\2 327r 



T Js~ 45 V'^M^I'^-^ ^^-^^ 
where Mp{= 1.2 x 10^^ GeV) denotes the Planck scale, and the right hand side is to be 
evaluated when the scale k~^, corresponding to the present horizon size, crossed outside the 
horizon during infiation. Equation (4.3) can be rewritten as 
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\^T/T\s ~ 0.067\/AAr^/^| ln(0^/M2)|i (4.4) 

where Nh{^ 55) denotes the number of e foldings experienced by this scale, and (pn is the 
field value when the scale crossed outside the horizon. 

The logarithmic factor in (4.4) is of order 1 - 10, assuming that the vacuum energy density 
that drives inflation is comparable to M^. For Mx ~ 10^^-^ GeV, M Mp = 1.2x 10^^ GeV, 
and taking (AT/T)cobe ~ 6 x 10^^, the fundamental quantity A is estimated to be 

A ^ 2.2 X 10^^^ (4.5) 

Any inflationary scenario is incomplete without an explanation of the origin of the ob- 
served baryon asymmetry in the universe. In the present case the inflaton mass ~ 
10^^'^ GeV, and so the basic idea is to create an initial lepton asymmetry via inflaton decay 
into one or more species of the heavy ('right handed') majorana neutrinos. The appearance 
of 'sphaleron' induced processes at the electroweak scale converts a specified fraction of this 
asymmetry into the observed baryon asymmetry. Details of how a satisfactory scenario is 



realized along these lines can be found in ref. |61]. We mention here a few salient features: 



(i) The reheat temperature ~ lO*^'^ GeV so that the out of equilibrium 
condition on the 'heavy' (~ 10^^ GeV) neutrinos is readily satisfied; 

(ii) The requirement that for temperature below the rates for lepton number 
violating processes {uu ffo*j^o* ^jjo ^ where H° is the electroweak 
scalar higgs) be smaller than the expansion rate {H ^ 20T^/Mp, where H 
denotes the Hubble constant) of the Universe, leads to the following constraint 
on the light neutrino masses ||6^ : 



4 eV 

< ^ 20 eV (4.6) 

(r,/10io GeV)5 



This is fully consistent with a cold plus hot dark matter scenario where neutrinos 
in the mass range 3 — 10 eV are needed. 
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(iii) In the present approach the colored scalar triplets which mediate proton 
decay are not needed for baryogenesis and consequently are allowed to have 
masses ~ Mx- 

(iv) Finally, it is possible, following ref. [Q, to identify the inflaton with the 
field that spontaneously breaks the axion symmetry. This would make for a more 
economical approach. 

(v) With an appropriate re-interpretation, the chaotic inflationary scenario can 
be realized within the framework of this 5*0(10) model. The ratio of the scalar 
to the tensor contribution in this case is 

(AT/r)|/(AT/r)| ^ 0.22. (4.7) 

B. Supersymmetric Inflation 

The presently measured gauge couplings of the standard model, when extrapolated to 
higher energies with supersymmetry (SUSY) broken at scales around 10'^ GeV 0], appear to 
merge at scales around 10^^ GeV . This is a boost for supersymmetric GUTS, with 5'f/(5) or 
50 (10) being the obvious gauge groups. In the presence of unbroken matter parity, either of 
them can provide a cold dark matter candidate in the form of LSP (lightest supersymmetric 
particle). However, hot dark matter in the form of massive neutrinos most naturally appear 
in the 5*0(10) model. The supersymmetric 50(10) scheme has some additional features 
which make it attractive from the particle physics viewpoint. For instance, 

(a) A Z2 subgroup of the center of 5*0(10) (more precisely Spin (10)) is left 
unbroken if tensor representations are employed to do the symmetry breaking. 
This Z2 symmetry which is not contained in SU (3)c x SU {2)^x11 (l)y, acts 
precisely as matter parity! 

(b) In some versions of SUSY 5*0(10), the important parameter tan/?(= (f)^/(f)'^, 
the ratio of the two vevs which provide masses to 'up' type and 'down' type 
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quarks) is predicted to lie close to rrit/mb [|6^. One consequence of this is the 



identification of the 'bino' (the supersymmetric partner of the f/(l)y gauge bo- 



son) as the LSP, with mass ~ 200 — 300 GeV. 



(c) Fermion Mass Ansatzes have recently attracted a fair amount of attention 



and are most simply realized within the framework of SO(IO) 



To summarize, particle physics considerations as well as observations of large scale struc- 
ture which favor a cold plus hot dark matter scenario, together suggest SUSY 5*0(10) as an 
attractive way to proceed. Infiation, either 'new' or 'chaotic', can be implemented by intro- 
ducing a suitable singlet superfield. Remarkably enough, singlets are typically employed to 
achieve the breaking of the GUT symmetry, (without breaking SUSY) and we exploit one 
of them to induce infiation! 

Let $ denote the 50(10) singlet (infiaton) superfield, xix) are the higgs superfields in 
the 126(126) representations whose vevs provide Majorana masses to the right handed neu- 
trino, and lQi{i = 1,2,3) are the matter superfields. To simplify the discussion, we restrict 
attention to the sector involving an interplay only between these superfields. This allows us 
to discuss the salient features of the (chaotic) infiationary scenario including baryogenesis. 
Consider the renormalizable superpotential 



Note that $ $ under the matter parity contained in 5*0(10) (similarly X:X ^ X: Xi and 
16, ^ -16,). 



The superpotential W gives rise to a supersymmetric ground state in which (vevs refer 
to the scalar components of the superfields) 



W = a^{xx- Ml) +f$2^f$3 



(4.8) 



+7ijl6il6jX 



< > 



0, < 16i >= 



(4.9) 



<X> 



<x>Vo 
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It is clear that matter parity is unbroken and we expect the LSP to contribute to the cold 
dark matter component. 

Even though B — L is now broken at Mx ~ 10^^ GeV, the right handed 'tau' neutrino 
mass must be of order 10^^ — 10^^ GeV, if the 'light' tau is to be the dark matter component. 
The inflaton must be at least twice as heavy, and one simple way to have ~ 10^^ GeV 
is to arrange the coefficients a and M in (4.8) to be of order m^/Mx and respectively. 
The decay rate of the inflaton into right handed neutrinos is given by 



With ~ 10^^ GeV , the reheat temperature T.^ is of order 10^ GeV . Baryogenesis via 



low Tr, the otherwise vexing gravitino problem is neatly avoided in this approach. 

Depending on the details, the spectral index n lies between 0.94 (if the quartic potential 
dominates) and 0.96 (with a quadratic potential dominant during the chaotic inflationary 
phase). The ratio (AT/T)^/ {AT/T)l ~ 0.22(0.11), respectively. 



The question we wish to address here is the following: Is it possible to implement inflation 
with GUTS without the gauge singlet? Surprisingly perhaps ||68|, an affirmative answer 
appears possible for a special class of supersymmetric GUTS in which, up to a normalization 
constant, the GUT scale is determined in terms of Ms and Mp, the SUSY breaking scale 
and the Planck scale respectively. Moreover, the normalization constant is fixed from the 
quadrupole anisotropy. Such models ||69| naturally arise after compactification of the ten 
dimensional Eg x Eg heterotic string theory ||70[|, and models based on G = SU{3)c x 
SU{3)l X SU{3)r or its subgroups provide some elegant examples. The scalar fields needed 
to spontaneously break G to SU{3)c x SU{2)l x U{1) can be used to drive infiation! 

The key ideas are relatively straightforward and perhaps best illustrated by a simplified 
example. Consider a rank five gauge symmetry H = SU{3)c x SU{2)i x f/(l)y x t/(l)', 




(4.10) 



leptogenesis now proceeds along the lines given in ref. Note that because of the relatively 



C. Inflation Without the Singlet 
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where the extra factor C/(l)' is to break at some superheavy scale M (below Mp). Let 
0, 4> denote the pair of higgs scalars whose vevs do this breaking. Note that < > = 
< >* so that the D term vanishes. Since the only independent dimensionful parameters 
are Ms and Mp (M is determined in terms of them), the terms in the superpotential 
are either absent or carry coefficients of order Ms- (Here $ denote the corresponding 
superfields.) Moreover, in order to ensure F flatness, the cubic couplings are also 

absent (otherwise 17(1)' would break at scales ~ Ms). 

In the absence of SUSY breaking the superpotential W is taken to be 

W^hxx^ + ^mf + --- (4-11) 

Ivi p 

where x denotes some matter superfield with the coupling h of order unity. Assuming a 
radiative breaking scenario along the lines envisaged in supergravity models with electroweak 

breaking, the effective potential takes the generic form 

y(0,0)~-M||0r + «:^^ (4.12) 
Minimization of (4.12) leads to the result 

M=|<0>| = |<0>|~ K~^{MsMp)^ (4.13) 

Provided that k!^ is sufficiently small, the potential in (4.12) will yield a satisfactory 
(chaotic) inflationary scenario. It turns out that k, ~ 10~^ (from COBE), which gives 
M ~ 10^^ GeV. The spectral index n ^ 0.92, while the ratio of the tensor to the scalar 
quadrupole anisotropy is (AT/T)|./(AT/T)| pa 0.4. Values of n ~ 0.88 (but no lower!) can 
be entertained within this framework. This is just as well since our analysis in the earlier 
sections seems to favor the range 1.0 > n > 0.9. 

To summarize, grand unification provides an elegant framework for implementing both 
the 'new' and 'chaotic' inflationary scenarios. Some popular models based on 5*0(10) or 
SU{3)c X SU{3)l X SU{3)r predict the value of the spectral index n in the range 0.96 
to 0.92. Cold dark matter, in axions and/or LSP, as well as hot dark matter in massive 
neutrinos are readily incorporated in these schemes. 
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V. CONCLUSIONS 



We have performed a ^ goodness of fit test of the predicted linear theory power spectra 
against data on scales ranging from 1 to 10^ Mpc, mainly from the COBE satellite and the 
QDOT IRAS survey of galaxies. We find that the inflation based scenario of large scale 
structure formation, in which the dark matter consists of cold plus hot components, can 
provide a good fit to large scale structure data. 

Taking the primordial power spectrum to have spectral exponent n, we find with 95% 
confidence, respectively, that (n^O.85), in models with (without) significant gravity 

wave contributions to the COBE anisotropy. The precise bound depends slightly on the 
HDM fraction. We find in models with a significant gravity wave anisotropy, that the 
COBE signal must not be dominated by the gravity wave contribution. 

If one insists on only one type of dark matter, i. e., CDM, the best fits are for n — 0.84 
[n = 0.92) in models without (with) gravity waves, respectively. 

The best fit region for all data, including some constraints from non-linear structure, is 
an roughly an ellipse (see figures 8 and 9). For models with small amplitude gravity wave 
anisotropics, the focii of the ellipse are approximately at {ytHDM = 0.20, n = 0.92) and 
{^HDM — 0.35, n — 0.98). For models with a large tensor COBE anisotropy, the ellipse 
is more eccentric, with the focii at roughly {^hdm = 0.2, n = 0.96) and {^hdm = 0.35, 
n = 0.98). Thus, the best fits occur for VLhdm ~ 0.20 — 0.35, with n very close to unity. 

Realistic examples of infiation from grand unification theories, including both supersym- 
metric and ordinary GUTS, which have these properties have been presented. These models 
are also consistent with other cosmological and particle physics constraints. 
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APPENDIX A: TRANSFER FUNCTION CALCULATIONS 

As described in ref. 0, we integrate the Fourier space evolution equations in (conformal) 
time using the gauge-invariant variables for the density Aca and velocity Va perturbations 
in each energy density component (CDM, neutrinos, photons, baryons) as given in ref. []7T 



Here we integrate the equations in conformal time instead of the scale factor as we had done 
previously. We begin the integration well before the matter dominated era {z = 8.3 x 10^ 
) and integrate up until the present time {z = 0). We have used only one flavor of massive 
neutrino, with the other two flavors essentially massless. For the baryon equations we 
use the equations for the difference between the baryon and photon density and velocity 
perturbations (the variables Sbr = ^c,baryon — (3/4)Acr and Vbr- = Vb — Vr given in ref. |7ll 



(section II-5). For the massive neutrinos we use the imperfect fluid treatment of ref. J72 

We numerically integrate the equations using a Haming predictor-corrector routine as we 

have found it tracks the oscillations of the relativistic components more accurately than 

Bulirsch-Stoer or Runge-Kutta routines. After recombination is completed, {z = 900) we 

switch from integrating the baryon-photon difference equations to simply integrating the 

baryon and photon component {Ac^baryon, ^cr and Vbaryon, Vr) equations separately. 

We have checked the results of our code by comparing our baryonic transfer functions 

against those of ref. @] who gives values for Qbaryon = 0.1 and 0.01 and found good agreement, 

although here we present results only for Qbaryon = 0.05. We have flt the transfer functions 

to an inverse 5th order polynomial and the results are given in table 2. We flnd a 5th 

order inverse polynomial works a little better for C+HDM models than a fourth order as 

is more usual for CDM models. The transfer functions given here are not baryonic transfer 

functions, but rather are flts to the total density perturbation {i.e., A = VLcdmAc^cdm + 

VLyAcu + ^baryonAc,baryon- Wc dcfluc our transfer function as 

A(fc,to)A(fc = 0,tO 

A{k,U)Aik = 0,to) ^ ' 
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where ti and to are the initial and present times, respectively. The transfer functions are 
accurate to a few percent down to A; = 1 /i/Mpc. 

^^^^ " 1 + tiA;0-5 + t2k^ + tgfci-s + Uk^ + hk^-^ ' ^^^^ 
In table 2, all coefficients are for k in Mpc (with h=0.5). 



APPENDIX B: ESTIMATING THE CONTRIBUTION OF NON-LINEAR 

EFFECTS. 

In the previous section we alluded to the fact that non-linear effects were becoming 
important at 8/i~^ Mpc. We would like to estimate the size of the non-linear effects. To 
get a crude estimate we use the "spherical collapse model" treatment (see, e.g., ref. [(73(1). 
This approximates a spherical overdensity in a flat universe locally as a miniature closed 
collapsing universe, so one can follow the collapse into the non-linear regime. We can define 
the non-linear overdensity as cr^^o^^i™ which is given by the following equation 

non-lin P . 9 (g - slug)^ 

= ;^-^ = 2(l-cos^)3' 
where ph is the background density and 9 is the conformal time coordinate which 
parametrizes a closed spacetime. For the same value of the linear theory predicts o"'-^™^ 

a^''^'^ = ^Hd - smd)f\ (B2) 

Thus if a(non-iin) = 0.6, we can estimate cx^'"^ = 0.4, a 50% correction. To keep within a 
range where the perturbations are linear we must restrict ourselves to scales where a < 0.4, 
where non-linear corrections are estimated to be < 30%. 

To estimate what value of am the counts in cells at £ = 20/i~^ Mpc implies, we must first 
correct for redshift space effects 



2 1 
367^ 56f 



al{l) (B3) 



We will assume for now that hi = 1.2, which implies am = 0.4. Thus we should not consider 
data on scales <20h~^ Mpc. 
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FIGURES 

FIG. 1. Transfer functions for the cold plus hot dark matter models. The curves represent 

the relative growth of density fluctuations as a function of scale. The curves represent the present 
time (z=0) transfer functions and in a universe with flcDM + ^hdm + ^baryan = !> where we have 
taken the canonical value of ^baryan = 0.05 from primordial nucleosynthesis. 

FIG. 2. We show the calculated model power spectra using our best fit model parameters for 
n = 1.00 only. The values of bg and bj have been fit as described in the text. The values of 6/ are 
close to 1.0, with the values bj = 1.1, 1.1, 0.9, 0.9, and 1.0 for nnoM = 0.0, 0.1, 0.2, 0.3, and 0.4, 
respectively. 

FIG. 3. The counts in redshift space cells data from the IRAS QDOT survey [25]. Illustrated 
are the curves for some selected models. Low mass fluctuation amplitudes of some models are 
compensated by using high values of the IRAS bias bj as can be seen here. 

FIG. 4. The values of the bulk streaming velocities extracted by the POTENT analysis 
program with 1 a error bars [29]. Also shown are predicted velocity curves for selected models. 
The models with low n and gravity waves can be ruled out by this data even though the cosmic 
variance for the velocity predictions in a single region are quite large. To illustrate this we have 
plotted the 68% confidence limits on the prediction of large scale streaming velocities for the 25% 
HDM model. 

FIG. 5. Here we show the IRAS QDOT power spectrum data [27] and the COBE amplitude 
constraint [23] (converted to a redshift space power spectrum constraint for b[ = 1.0). The mod- 
els shown have negligible gravity wave contributions to the COBE anisotropy, and have best fit 
normalizations bs = 1.26, 1.41, and 1.54 for the models shown with 0%, 25%, and 45% HDM, 
respectively. We have shown only half of the QDOT data so the plot is easier to read. The value 
of n = .84 yields the best fitting CDM model without gravity waves. Note that simply tilting a 
CDM spectrum cannot reproduce the large scale "bump" in the IRAS power spectrum, which is 
generally why CDM models do not fare as well as C+HDM models. 
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FIG. 6. Contour plot for - No gravity waves contours in the n-Q,HDM plane for models 
with a negligible gravity wave content. Moving outward from the center of the graph are contours 
corresponding to .5 % (dotted), 1% (dash-dotted), 5% (dashed), 10% (heavy dashed) 25% (long 
dashed), 50% (heavy long dashed), 68% (solid), and 95% (heavy solid) confidence levels. 

FIG. 7. Contour plot for x^ - With gravity waves x^ contours in the u-Uhdm plane for models 
with gravity waves. Moving outward from the center of the graph are contours corresponding to .5 
% (dotted), 1% (dash-dotted), 5% (dashed), 10% (heavy dashed) 25% (long dashed), 50% (heavy 
long dashed), 68% (solid), 95% (solid), 99% (heavy solid) confidence levels. 

FIG. 8. Same as figure 6, but we have added the constraints from non-linear data. We have 
added the restriction that bg > 1.25 and that we have sufficient power for quasar formation. The 
procedure for adding these constraints in a x^ analysis is described in the text. The significance of 
the contours is here not well defined because of the way we have added the non-linear constraints. 

FIG. 9. Same as figure 7, but we have added the constraints from non-linear data. We have 
added the restriction that > 1.25 and that we have sufficient power for quasar formation. The 

significance of the contours is here not precise because of the way we have added the non-linear 
constraints. 
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TABLES 



TABLE L /.', f f i:)r(il)od by couuls in colls o[ A"olumo f'' 



£ keff 


20 Mpc/h 


0.120 h/Mpc 


30 Mpc/h 


0.081 h/Mpc 


40 Mpc//i 


0.031 /i/Mpc 


60 Mpc/h 


0.022 h/Mpc 



TABLE II. transfer functions with ^baryon = 0.05, h = 0.5 





h 


t2 


t3 


t4 


h 


0.00 


-1.150 


29.60 


48.49 


-43.17 


132.4 


0.05 


-0.8654 


17.65 


165.1 


-277.1 


343.4 


0.10 


-0.2942 


1.393 


274.6 


-472.6 


538.2 


0.15 


0.1157 


-8.820 


330.0 


-541.4 


660.8 


0.20 


0.3176 


-12.69 


334.9 


-495.8 


726.1 


0.25 


0.3128 


-10.60 


296.3 


-361.6 


756.5 


0.30 


0.1363 


-3.540 


219.1 


-142.9 


771.2 


0.35 


-0.1454 


6.144 


127.6 


78.64 


840.3 


0.40 


-0.4276 


15.15 


53.35 


193.1 


1084. 


0.45 


-0.6522 


21.44 


15.98 


131.2 


1582. 


0.50 


-0.7882 


24.01 


25.49 


-117..'3 


2395. 
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